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ABSTRACT 

Near-infrared spectra of four QSOs located at z > 6 are obtained with the OH- 
airglow suppressor mounted on the Subaru telescope. The Fe Il/Mg II emission- line 
ratios of these QSOs are examined by the same fitting algorithm as in our previous 
study of z < 5.3 QSOs. The fitting results show that two out of the four z > 6 QSOs 
have significant Fe II emission in their rest-UV spectra, while the other two have almost 
no Fe II features. We also applied our fitting algorithm to more than 10,000 SDSS QSOs 
and found two trends in the distribution of Fe Il/Mg II against redshift: (1) the upper 
envelope of the Fe Il/Mg II distribution at z > 3 shows a probable declination toward 
high redshift, and (2) the median distribution settles into lower ratios at z ~ 1.5 with 
small scatter compared to the other redshift. We discuss an Fe/Mg abundance evolution 
of QSOs with a substantial contribution from the diverse nature of the broad-line regions 
in high-redshift QSOs. 

Subject headings: galaxies: active — quasars: emission lines — quasars: general — infrared: 
general 



1. INTRODUCTION 

The Fe/Mg abundance ratio, one of the 
most important parameters for estimating the 
age of a stellar system after its initial star- 
burst, is based on the delay of iron enrich- 
ment by Type la supernovae and compared 
with the formation of such a-elements as mag- 
nesium by Type II supernovae (see Hamann 
& Ferland 1999 for a review). The delay of 
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iron formation was generally thought to be 
~ 1 Gyr (e.g., Yoshii et al. 1996), while more 
recent studies suggest that this delay can be 
~ 0.3 Gyr, depending on environmental con- 
ditions of the star formation process (Friaca 
& Terlevich 1998; Matteucci & Recchi 2001). 
The Fe Il(UV)/Mg II emission-line ratio of 
QSOs is considered a probable indicator of the 
Fe/Mg abundance ratio, which has been mea- 
sured in various redshift ranges: 3.1 < z < 4.7 
by Thompson et al. (1999, hereafter THE99), 
< z < 5.3 by Iwamuro et al. (2002, here- 
after Paper I), 5.7 < z < 6.3 by Freudling 
et al. (2003), < z < 4.8 by Dietrich et al. 
(2003), and 3.0 < z < 6.4 by Maiolino et al. 
(2003). The median values of the Fe Il/Mg II 
emission-line ratios are almost constant at all 



redshifts, while large scatter dominates the 
distribution of the ratios at z > 3, indicating 
the diversity of the formation histories of high- 
redshift QSOs. On the other hand, recent 
photoionization calculations of Fe + atoms in 
the broad-line regions of QSOs show that the 
Fe Il(UV)/Mg II emission- line ratio is sensi- 
tive not only to Fe and Mg abundance but 
also to physical conditions, including micro- 
turbulence (Verner et al. 2003). However, the 
median value of the Fe Il/Mg II emission-line 
ratio is expected to converge on a significantly 
smaller value when the age of the universe be- 
comes less than ~ 0.3 - 1 Gyr, assuming that 
the substantial scatter caused by the diverse 
characteristics (or microturbulence velocities, 
etc.) of the broad-line regions in QSOs has 
the same contribution at any redshift. 

In this paper, we report the Fe II/Mg II 
emission-line ratios of four QSOs at z > 6 ob- 
served by the OH-airglow suppressor (OHS; 
Iwamuro et al. 2001) and CISCO (Motohara 
et al. 2002) mounted on the Subaru telescope. 
In §2 we analyze the data by the same method 
as in Paper I, and then we compare the re- 
sults with other observations in §3. In §4 we 
discuss the expected Fe/Mg abundance evo- 
lution of QSOs using the combined samples 
of these results, the new Sloan Digital Sky 
Survey (SDSS) archival data 1 , and the previ- 
ous samples in Paper I. Throughout the paper 
we adopt a cosmology with Hq = 70 km s _1 
Mpc- 1 , Q M = 0.3, and Q A = 0.7. 

2. OBSERVATIONS AND DATA RE- 
DUCTION 

The observations were carried out on 2002 
February 27 and March 2 and again on 2003 
March 20 and 21, using OHS and CISCO 
mounted on the Subaru telescope. Since the 
Fe II and Mg II are redshifted to the H and 
K bands, respectively, for z > 6 QSOs, the 

1 See http://www.sdss.org/drl. 



K-band spectra were obtained by CISCO and 
the Jii-band spectra by OHS separately in 
each observation run. The typical exposure 
sequence is 200 s x 3 x 4 positions (CISCO) 
or 1000 s x 4 positions (OHS), in which the 
object is moved about 10" along the slit by 
nodding the telescope after every three expo- 
sures (CISCO) or one exposure (OHS). The 
slit widths were 0"8 (CISCO) and 1" (OHS), 
providing spectral resolutions of 330 and 210, 
respectively. After this exposure sequence, a 
nearby SAO star with a spectral type of A or 
F was observed as a spectroscopic standard 
to remove the telluric atmosphere absorption 
features and to correct the instrumental re- 
sponse. The pixel scale was O'.'106/pixel with 
the infrared secondary mirror of the telescope, 
while the seeing size was 0'.'7. To calibrate 
the relative flux between the K- and Ji^-band 
spectra, short imaging exposures of 20 s x 3 x 
4 positions were executed by CISCO in both 
the K' (1.97-2.30 /an) and H (1.50-1.79 /an) 
bands, except for SDSS 1030+0524. The ob- 
servation log is summarized in Table 1. 

The obtained data were reduced using 
IRAF with typical reduction procedures for 
infrared images (see Paper I). The resultant 
one-dimensional spectra in the K and JH 
bands were combined using the photometric 
results of the imaging observation by CISCO. 
As a result, through a circular aperture of 2'/ 2 
diameter, the average fluxes between 1.97- 
2.30 and 1.50-1.79 /im correspond to the K'- 
and .ff-band magnitudes respectively. The ex- 
traction procedure for Fe II and Mg II emis- 
sion lines from the combined spectra is the 
same as in Paper I: nonlinear \ 2 fitting for the 
Mg II (Moffat function) and for the continuum 
(power law) profile with six free parameters, 
after subtraction of the Fe II template from 
the composite QSO spectrum of the Large 
Bright Quasar Survey (LBQS; Francis et al. 
1991) and the Balmer continuum template 
from the UV spectrum of 3C 273 (Wills et al. 
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1985) with various multiplying factors (see eq. 
(l)-(4) in Paper I). The fitting range is 2150- 
3300A, which is also the wavelength range for 
our definition of Fe II flux. Note that the 
strength of the Fe II emission is determined 
by shape comparison between our Fe II tem- 
plate and the observed spectral features (large 
bumps and small depressions), which does not 
require the assumption of the power-law con- 
tinuum level (we need the assumption of the 
shape of the Fe II template instead). 

Although the Moffat function is very useful 
for fitting the various profiles of Mg II emis- 
sion lines with the minimum number of pa- 
rameters, a serious problem is that the fitted 
Mg II profile sometimes has very broad wings. 
To prevent such unexpected broad wings, we 
set the lower limits at b > 0.5 for the Moffat 
power (eq. (1) in Paper I) and reanalyzed all 
the data fitted with b < 0.5 in Paper I (see 
Appendix). 

3. RESULTS 

The object spectra with fitted components 
are shown in Figure 1, and the numerical 
results are listed in Table 2. The spectra 
with wavelengths longer than 1.80 /im are ob- 
tained by CISCO, and the shorter parts by 
OHS. The strength of the Fe II emission is 
mainly determined by the break feature at 
rest 2200 A, which is observed by OHS with su- 
perior signal-to-noise ratios. The wavelength 
range between 1.80 and 1.95 /im corresponds 
to the atmospheric absorption band, whose 
contribution to the fitting results is very small 
because of large errors. 

As shown in Figure 1 and Table 2, we de- 
tected strong Fe II emission from SDSS J1148+5251, 
which is consistent with Maiolino et al. (2003). 
SDSS J1048+4637 also shows significant Fe II 
emission, while the Fe Il/Mg II ratio is not 
as much as the value of ~8.1 reported by 
Maiolino et al. (2003). The absorption-line 



feature at 1.29 /im in SDSS J1048+4637 
is the broad absorption line component of 
C III] reported by Maiolino et al. (2004). 
The observed spectra of these objects fit well 
into our template spectrum of Fe II in ev- 
ery detail. Although Barth et al. (2003) 
reported a smaller Fe Il/Mg II ratio of 
~4.7 for SDSS J1148+5251, this inconsis- 
tency is mostly caused by the difference in 
the Fe II template under Mg II emission. 
On the other hand, SDSS J1030+0524 and 
SDSS J1630+4012 show almost no Fe II emis- 
sion or any break feature at rest 2200A. These 
results for SDSS J1030+0524 are consistent 
with Freudling et al.(2003), as well as the un- 
known absorption-line feature at 1.57 fim. As 
a result, even at z > 6, the Fe Il(UV)/Mg II 
emission-line ratios still show large diversity. 

4. DISCUSSION 

Figure 2 shows the Fe Il/Mg II ratio and 
the absolute AB magnitude at rest 2500A as 
a function of redshift, for which we reanalyzed 
the SDSS data on the basis of the DR1 quasar 
catalog (Schneider et al. 2003). The other 
z < 5.3 data are the same as those plotted 
in Figure 7 of Paper I, except for four objects 
listed in the Appendix. The median values 
and the standard deviation of the sample dis- 
tribution are shown by the squares with the 
cross bars in this figure whose numerical val- 
ues are listed in Table 3. Although the num- 
ber of samples at z > 3 is insufficient to make 
a statistical argument, we can deduce the 
trends of the evolution of the Fe Il/Mg II ratio 
from Figure 2. First, the median Fe Il/Mg II 
ratios at z > 3 are almost constant, while the 
upper envelope of the distribution plotted by 
the dashed line in Figure 2 shows a probable 
declination toward high redshift. Second, the 
median distribution settles into lower ratios at 
z ~ 1.5, with small scatter compared to the 
other redshift. Here, the downward error bars 
from these median points are almost equal to 
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the typical fitting errors of the samples in each 
bin (column 3 in Table 3) , while the longer up- 
ward error bars are affected by the substantial 
scatter of the Fe Il/Mg II ratios. 

Assuming that the Fe/Mg abundance ra- 
tio affects these rough trends, these small 
Fe II/Mg II ratios at z ~ 1.5 originated with 
the dilution of the Fe abundance with the out- 
flow gas from low-mass stars after ~3 Gyr 
from the initial starburst (Yoshii et al. 1998). 
The scatter of the sample at this redshift (a 
factor of ~2 between the top and bottom 
of the cross bars) corresponds to the max- 
imum contribution of the scatter caused by 
the differences in the physical conditions of 
the broad-line regions in high-redshift QSOs 
(Verner et al. 2003). If this substantial scatter 
of factor of ~2 is the universal nature of the 
QSOs, independent of the redshift, the larger 
scatter of the z > 3 samples are affected by 
the difference in the Fe/Mg abundance ratio 
or the time passage after the initial starburst. 
This idea is also consistent with the declina- 
tion of the upper envelope of the Fe II/Mg II 
distribution toward high redshift, which is ex- 
pected to be the initial abundance evolution 
of QSOs. 

An alternative explanation for the varia- 
tion of the Fe II/Mg II ratios at z ~ 2 may 
be their luminosity dependence such that the 
QSOs of higher luminosities, in a low-redshift 
(z =0.1-0.6) sample from Hubble Space Tele- 
scope (HST) UV archives (Tsuzuki 2004) 
as well as in a high-redshift (z =3-5) sam- 
ple from near-infrared observations (Iwamuro 
et al. 2002; Dietrich et al. 2003), are found to 
have higher Fe II/Mg II ratios. The larger 
variation of the Fe II/Mg II ratios for the 
brighter QSOs suggests that either the ion- 
izing photon flux arriving at the broad-line- 
emitting clouds or the physical characteris- 
tics of the clouds would scatter more signif- 
icantly for brighter QSOs. To estimate their 
respective contributions to the Fe II/Mg II ra- 



tio is a theoretical challenge that requires an 
elaborate photoionization model of Fe + in the 
broad-line regions of QSOs. 
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A. Re-analyzed Objects Having Broad Mg II Wings 

We re-analyzed all the data having broad Mg II wings (b < 0.5) in Paper I with lower limits at 
b > 0.5 for the Moffat profile of 



FMgllW = i 



1 + 



A 



27981 V 
a J 



(Al) 



All the data were fitted with b = 0.5 (Lorentzian), and the results are listed in Table 4. Three out 
of seven objects listed in Table 4 are not plotted in Figure 2 because we adopted the results of the 
OHS data rather than those of THE99 and the results of the HST FOS archive rather than those 
of Kinney et al. (1991) (see Table 2 and Table 4 in Paper I). 
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Fig. 1. — Observed spectra of the four QSOs listed in Table 1. A power law continuum, Mg II, Fe II, 
a Balmer continuum (including other lines except for Mg II), and the sum of these components are 
plotted for comparison with the observed spectrum with the boxcar-smoothing of six pixels. The 
two vertical lines indicate the fitting range of 2150-3300A in the rest-wavelength. 
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Fig. 2. — Redshifts vs. Fe Il/Mg II and absolute magnitudes. The fitting results for the QSO 
spectra from Kinney et al. (1991), H ST /FOS and SDSS archives (dots), THE99 (triangles), our 
data in Paper I (open circles), and this work (filled circles) are plotted with the median values 
(squares) of the appropriate bins. The horizontal solid line indicates the Fe Il/Mg II of 4.09 for 
the LBQS composite spectrum (see Paper I). The squares with the cross-bars represent the median 
values and the standard deviations of the sample distribution, whose upper envelope is connected 
by a broken line. We assume Hq = 70 km s™ 1 Mpc" 1 , = 0.3, and = 0.7 for the estimation 
of M2500 (the absolute AB-magnitude at rest-frame 2500 A) and the cosmic age at the redshift. 
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Table 1 
Log of Observations 



Exposure Coordinate 



Object Name 


Redshift 


if'-mag a 


i7-mag a 


Date 


Instrument 13 


Time (s) 


Seeing 


Reference 


SDSS J103027.10+052455.0 


6.28 


17.67 


c 


2002 Feb. 27 


CISCO 


2400 


0780 


1 








18.57 


2002 Mar. 2 


OHS 


4000 


Of.' 88 




SDSS J104845.05+463718.3 


6.23 


17.12 


17.76 


2003 Mar.20 


CISCO 


2400 


0756 


2 








17.83 


2003 Mar.21 


OHS 


4000 


Of.' 71 




SDSS J114816.64+525150.3 


6.43 


16.98 


17.70 


2003 Mar.20 


CISCO 


2400 


Of.' 65 


2 








17.62 


2003 Mar.21 


OHS 


4000 


0762 




SDSS J163033.90+401209.6 


6.05 


18.40 


19.25 


2003 Mar.20 


CISCO 


2400 


0771 


2 








19.18 


2003 Mar.21 


OHS 


8000 


0768 





a Observed magnitudes by CISCO and OHS short-imaging observations with a 272 circular diameter aperture. Typical photometric 
errors are 0.05 mag. 

b The slit widths are 078 (CISCO) and 1" (OHS), corresponding to spectral resolutions of 330 and 210, respectively. 
c The -ff-band imaging observation was not executed. 



References. — (1) Fan et al. (2001); (2) Fan et al. (2003). 



Table 2 

Results of Fitting Calculations 



FWHM a /(FeII) b J(Mg II) EW(Fe II) b £W(Mg II) 



Object Name z (km s" 1 ) (10~ 16 ergs s" 1 ™" 2 ) Fe II/Mg II b (A) (A) 



SDSS J103027.10+052455.0 6.311 3590 S^l^g 6.04±1.40 0.99^^ 21.0 22.9 

SDSS J104845.05+463718.3 6.203 4050 42.4±14.7 11.3±2.1 3.74±l'.47 81.8 25.7 

SDSS J114816.64+525150.3 6.403 3020 137±16 14.4±2.2 9.52±1.82 333 41.7 

SDSS J163033.90+401209.6 6.065 2690 7.30±7.00 7.12±1.61 1.02±1.01 52.2 55.8 



a FWHM of the Mg II emission-line in the rest-wavelength. 
b Fe II is defined over the domain of 2150-3300A. 



Table 3 
Median values of Fe II/Mg II 



Redshift a Fe II/Mg II b Error c Number d 



0.101±0.068 4.48l£jjj| 

0.839±0.089 3.36±J;|| 

1.003±0.075 3.49±g;|3 

1.144±0.066 3.34lg;§| 

1.273±0.063 3.05±g;|| 

1.402±0.066 2.78lJ;§? 

1.528±0.060 2.79l5;fg 

1.647±0.060 2.81±5;§| 

1.769±0.063 3.00l5;|? 

1.907±0.075 3.75li^ 

2.091±0.109 5.11±f;gg 

3.411±0.220 5.50l|^ 

4.475±0.215 4.40±J§? 

4.997±0.291 5.69±|:?f 

6.234±0.169 2.381^7 



0. 


.84 


30 


0. 


.82 


1101 


0. 


90 


1101 


0. 


.90 


1101 


0. 


.92 


1101 


0. 


93 


1101 


0. 


.81 


1101 


0. 


.84 


1101 


0. 


.95 


1101 


1. 


.24 


1101 


1. 


.67 


1099 


1. 


.45 


6 


1. 


.76 


9 


2. 


.70 


8 


2. 


.02 


4 



a Redshift range of each bin. 

b Median value with standard deviation of the 
sample distribution. 

c Typical (root mean square) fitting errors of 
Fe II/Mg II ratios of the samples included in each 
bin. 

d Number of samples. 
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Table 4 

Results of Fitting Calculations for THE99 data 







FWHM 




EW(Fe II) 


EW(Mg II) 




Object Name 


z 


(km s -1 ) 


Fe II/Mg II 


(A) 


(A) 


Reference 


BR 1033-0327 


4.527 


4050 


2.34±0.78 a 


98.5 


46.6 


THE99 


1358+391 


3.288 


4610 


2.86+0.49 


110 


46.1 


THE99 


PKS 2126-158 


3.290 


3430 


2.92+0.63 


71.3 


29.1 


THE99 


IIIZw 2 


0.089 


3020 


1.95±0.28 a 


155 


79.1 


IUE 


PG 0844+349 


0.065 


2320 


6.31+0.23 


238 


40.4 


HST 


IRAS 1334+24 


0.109 


3010 


2.62+0.27 


138 


40.9 


HST 


MRC 2251-178 


0.064 


4200 


3.53±0.27 a 


424 


143 


IUE 



a These data are not plotted in Figure 2 (see Appendix). 

References. — (IUE) Kinney et al. (1991) (ftp://dbc.nao.ac.jp/DBC/NASAADC/catalogs/3/3157/); 
(HST) HST archive in CADC (http://cadcwww.dao.nrc.ca/hst/science.html); (THE99) 
Thompson et al. (1999) 
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